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ABSTRACT 

The transmission of genetic information relies on 
Watson-Crick base pairing between nucleoside 
phosphates and template bases in template-primer 
complexes. Enzyme-free primer extension is the 
purest form of the transmission process, without 
any chaperon-like effect of polymerases. This sim- 
ple form of copying of sequences is intimately linked 
to the origin of life and provides new opportuni- 
ties for reading genetic information. Here, we re- 
port the dissociation constants for complexes be- 
tween (deoxy)nucleotides and template-primer com- 
plexes, as determined by nuclear magnetic reso- 
nance and the inhibitory effect of unactivated nu- 
cleotides on enzyme-free primer extension. Depend- 
ing on the sequence context, IQ's range from 280 
mM for thymidine monophosphate binding to a ter- 
minal adenine of a hairpin to 2 mM for a deoxyguano- 
sine monophosphate binding in the interior of a se- 
quence with a neighboring strand. Combined with 
rate constants for the chemical step of extension 
and hydrolytic inactivation, our quantitative theory 
explains why some enzyme-free copying reactions 
are incomplete while others are not. For example, for 
GMP binding to ribonucleic acid, inhibition is a sig- 
nificant factor in low-yielding reactions, whereas for 
amino-terminal DNA hydrolysis of monomers is crit- 
ical. Our results thus provide a quantitative basis for 
enzyme-free copying. 

INTRODUCTION 

Step-wise extension of a growing oligonucleotide chain by 
nucleotides, directed by a template, is the molecular basis 
of replication and transcription (1,2). This process is found 
in all kingdoms of life. Its rate and fidelity are critical for 
the survival of species, both in terms of passing on genetic 
information, and in terms of allowing for mutations (3). 
Watson-Crick base pairing is known to be the guiding prin- 



ciple of nucleobase selection during replication, but numer- 
ous factors affect the rates of polymerase-catalyzed exten- 
sion (4-6). Some alternative base pairs are accepted by poly- 
merases (7-10), but it is not easy to dissect the contribu- 
tions that the active site and the template/primer duplex 
make to the interactions driving the incorporation of nu- 
cleotides. Some dNTP analogs are readily incorporated by 
polymerases, even though their base pairs destabilize du- 
plexes (11-13). 

One field where the strength of the template effect pro- 
vided by base pairing between nucleotides and templates 
is particularly important is enzyme-free replication (14). 
Enzyme-free or 'chemical' primer extension, is solely driven 
by the template effect experienced by an incoming nu- 
cleotide binding to primer-template complexes and the in- 
trinsic reactivity of the monomer (15). A quantitative un- 
derstanding of this reaction is critical for theories on how 
life might have arisen during an early phase of evolution 
(16,17). It is currently unclear whether the template effect 
provided by a template: primer complex is sufficient to al- 
low for enzyme-free copying of sequences long enough to 
act as polymerase ribozymes (18,19). 

Enzyme-free copying was first demonstrated for ribonu- 
cleic acid (RNA) (20-24). Detailed studies suggested that 
only sequences rich in cytidylic acid were able to induce the 
spontaneous formation of complementary strands (25), and 
the prospect of replication in systems containing all four nu- 
cleotides was called 'remote' (26). Later work showed that 
with oxyazabenzotriazole leaving groups, low temperatures 
and downstream-binding oligonucleotides acting as 'helper 
strands' rates are accelerated and yields improve (27,28). 
Further, a competitive inhibition by unactivated nucleotides 
(produced through hydrolysis in the time course of an assay) 
can be avoided when template and primer are immobilized 
and the supernatant containing the monomers is removed 
periodically (29). A similar approach has recently been im- 
plemented for vesicles, using periodic dialysis against solu- 
tions of activated monomers (30). In favorable cases, exten- 
sion by any of the four nucleotides (A/C/G/U) opposite 
their canonical base pairing partner occurs (29). This ap- 
proach requires intervention and well defined conditions, 
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Figure 1. Binding equilibrium between nucleotide and primer-template 
duplex; B, B' = nucleobases. 



though (31), and serious challenges remain for monomer- 
based self-replicating systems (32). 

Fast rates in enzyme-free primer extension can be 
achieved by using amino -terminal primers, whose exten- 
sion by reactions with activated monomers produces ge- 
netic polymers with phosphoramidate linkages that are iso- 
electronic to natural phosphodiesters (33-39). With amino- 
terminal primers, high yields were observed for any of the 
64 triplets whose center nucleotide provides the templating 
base (40). Further, slow extension of primers after a mis- 
incorporation has been shown to improve the fidelity of 
copying of longer sequences, thus avoiding a potential 'error 
catastrophe' in enzyme-free replication (41^3). What has 
remained unclear, though, is what fraction of the reactiv- 
ity observed in the extension of amino-terminal primers is 
due to template-independent chemical reactivity and what 
fraction is due to the template effect. Overextension of the 
primer, beyond the length of template, suggests that the con- 
tribution of non-templated polymerization can be signifi- 
cant (44). Without a quantitative understanding of the tem- 
plate effect, the question of whether longer sequences can 
be copied under prebiotically plausible conditions remains 
difficult to answer. A quantitative understanding of the tem- 
plate effect will also help to understand how much the sub- 
strates contribute to the fidelity of polymerase-catalyzed 
polymerization. Finally, such data will help to develop new 
methods for reading out genetic information in enzyme-free 
fashion (45,46). 

Base pairing between isolated nucleobases has been stud- 
ied in organic solvents (47,48). But, there appear to be no 
experimental binding constants for complexes between nu- 
cleotides and primer-template duplexes (Figure 1), even 
though there has been intense theoretical work (49-51). For 
example, Bickelhaupt and colleagues have calculated bind- 
ing energies for model complexes of template-primer du- 
plexes and incoming nucleoside phosphates. Their AG a ffi n i ty 
values range from —20.8 kcal/mol for an incoming G pair- 
ing with C as templating base, when intrinsic thermal and 
entropy effects are ignored, and —3.1 kcal/mol for an in- 
coming T pairing with A when a large estimated value for 
such effects is used. 

Here, we report binding constants for the complexes of 
nucleotides with primer-template duplexes, as determined 
by two complementary techniques. A quantitative model 
presented produces time-dependent yields of enzyme-free 
extension. Our data has allowed us, for the first time, to 
understand the concentration and sequence dependence of 



enzyme-free primer extension, as well as the inhibitory ef- 
fect of hydrolyzed monomers. 

MATERIALS AND METHODS 

Activated monomers and primers 

Oxyazabenzotriazolides (OAt esters) of nucleotides (dN- 
MPs or GMP) were prepared via activation of with 
EDC/HOAt (35). Primers with 3 / -terminal 3 / -amino-2 / ,3 / - 
dideoxynucleoside residue were prepared on as described 
previously (52) and were purified by HPLC. Unmodified 
DNA strands were purchased from Biomers (Ulm, Ger- 
many) in salt-free form and were used without further 
purification. Unmodified hairpin sequences and the RNA 
hairpin were purchased from Biospring (Frankfurt, Ger- 
many). Further details can be found in the Supplementary 
material. 

NMR experiments 

Nuclear magnetic resonance (NMR) samples (200 fxL) were 
prepared in 3 mm semimicro tubes and were 0.5 mM in the 
hairpin oligonucleotide. Spectra were recorded on a Bruker 
Avance 500 spectrometer. Signal assignment was based on 
a combination of NOESY and TOCSY spectra with pre- 
saturation to suppress the solvent signal. Data processing 
used an exponential function with a line broadening setting 
of 0.3 Hz. For representative two-dimensional (2D) spectra 
and more detailed protocols, please see the Supplementary 
material. 

Primer extension assays 

Primer extension assays with MALDI-ToF-based analy- 
sis were performed as previously described (40). For in- 
hibitor assays, experiments were typically performed as fol- 
lows. To the assay solution (10 |xL final volume) contain- 
ing the primer (36 |xM), the template (54 jxM) and, unless 
noted otherwise, the downstream-binding oligonucleotide 
(54 |jlM) in HEPES buffer (200 mM, 400 mM NaCl, 80 
mM, pH 8.9), the unactivated 3 / -deoxynucleotide (la-t) was 
added using an aliquot of a stock solution (72 mM) in as- 
say buffer. Assays were started by addition of an aliquot 
of the aqueous stock solution of the OAt-esters of 3'- 
deoxynucleotides (final concentration 3.6 mM for 7a, 7c or 
7g, and 7.2 mM for 7t). Additional details and data can be 
found in the Supplementary material. 

RESULTS 

Nucleotides binding to terminal template bases 

Two experimental systems were used to measure the bind- 
ing of nucleotides to templating bases (Figure 1). The first 
used NMR spectroscopy as monitoring technique and short 
hairpins as intramolecular models of primer-template du- 
plexes (Figure 2). The overhang at the 5 / -terminus of the 
hairpin provided the templating base to which deoxynu- 
cleoside monophosphates (dNMPs) la-t were allowed to 
bind. The stem of the hairpin oligonucleotides was chosen 
to be long enough to give a stable helix at room tempera- 
ture, but short enough to allow for NMR without isotope 
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Figure 2. Nucleotides and hairpins used for NMR titration. Loops are 
hexaethylene glycol linkers (HEG). 
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Figure 3. Typical results from NMR titrations, (a) Overlay of spectra with 
^-NMR signal of H-8 proton of the S'-terminal A residue of hairpin 5'- 
ATGC(HEG)GCA (2a) (0.5 mM) with increasing concentration of TMP 
(It) (0 to 774 mM) at 20°C. (b) Plot of chemical shift displacement upon ad- 
dition of TMP (It). See Supplementary Figures S8, S9 and S12-S14 (Sup- 
plementary material) for additional spectra and plots of chemical shifts. 



enrichment. A hexaethylene glycol (HEG) linker loop (53) 
was chosen to provide stability without complicating as- 
signment. 

A series of l H NMR spectra of hairpins 3a and 3c at 
increasing temperature in solutions with buffer conditions 
typical for primer extensions (35) confirmed that the melt- 
ing transition of the helices was above 40° C (see Supplemen- 
tary Figures S3-S5, Supplementary material). The assign- 
ment of the hairpin NMR signals was based on 2D spec- 
tra and a literature-known assignment strategy (54). Titrat- 
ing deoxynucleotides into the NMR solution of hairpins 
led to downfield or upheld shifts of resonances of the ter- 
minal residues of the hairpin (Figure 2 and Supplementary 
Figures S8-S9 in the Supplementary Material). Curve fit- 
tings to plots of shifts against nucleotide concentrations 
then gave binding constants (see Supplementary material, 
Supplementary Figures S12-S14 for details) (55,56) (Figure 
3). 
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Figure 4. Binding equilibria underlying the inhibitory effect of a free nu- 
cleotide on chemical primer extension. 



Binding constants for complementary pairs of nucleotide 
and templating base were found to range from 10 to 280 
mM, depending on the templating base, the length of the 
overhang, and the base at the 3 '-terminus of the primer seg- 
ment of the hairpin (Table 1). Large changes in were 
observed when switching from one templating base to an- 
other. The binding of thymidine 5 / -monophosphate (TMP, 
It) to adenine-displaying hairpins 2a, 3a and 4a is ~20-fold 
weaker than that that of dGMP (lg) to hairpins 2c and 
3c. The values for 2 / -deoxycytidine 5 / -monophosphate 
(dCMP, lc) and dAMP (la) binding to the hairpins display- 
ing their complementary base as templating residue (2g and 
4t) are in between those measured for TMP and GMP, with 
values of 38 and 40 mM, respectively. The base at the ter- 
minus of the hairpin helix has a modest effect on the 
A mismatch between incoming nucleotide and templating 
base leads to ~ 10-fold drop in K& for a G:T wobble pair 
(lg binding to 4t) and a more than 100-fold drop for a C:T 
base combination (lc and 4t). An exploratory measurement 
with an all-RNA system (5g:6c) gave a of 14 mM, which 
is close to that measured for the same G:C pairing in the 
DNA hairpins (lg:2c and lg:3c). 



Binding in the interior of longer sequences 

We then developed an approach for measuring binding 
of deoxynucleotides to the primer extension site of longer 
DNA templates in the presence or absence of downstream- 
binding oligodeoxynucleotides. We used the inhibitory ef- 
fect of free nucleotides added to primer extension mixtures 
for our measurements. The free nucleotide competes with 
the activated nucleotide for the primer extension site, thus 
inhibiting the reaction (Figure 4). 

Monoexponential fits to data sets from a series of kinet- 
ics runs with increasing concentration of unactivated de- 
oxynucleotide gave inhibitory constants (K\ n h), from which 
dissociation constants (^d) were calculated. The kinetic 
analysis was based on the following model. (A discus- 
sion of its implications and a more sophisticated model 
can be found in the Supplementary material.) The primer 
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Table 1. Dissociation constants for nucleotides binding to hairpin termini, as determined by NMR titration* 



Nucleotide 


Templating base 


Hairpin 




K d (mM) b 


dTMP (It) 


A 


3 / -ACG(HEG)CGTA 
Y-ACGiTTFG^CGTATC 


(2a) 


260 


dTMP (m\ 






280 


dTMP (It) 


A 


3 / -GTC(HEG)GACA 


(4a) 


240 


dGMP (lg) 


C 


3 / -ACG(HEG)CGTC 

3 / -ACG(HEG)CGTCTC 

3 / -GTC(HEG)GACT 


(2c) 


10 


dGMP (lg) 


C 


(3c) 


16 


dGMP (lg) 


T (wobble pair) 


(4t) 


100 


dCMP (lc) 


G 


3 / -ACG(HEG)CGTG 


(2g) 


40 


dCMP (lc) 


T (mismatch) 


3'-GTC(HEG)GACT 


(4t) 


>2000 d 


dAMP (la) 


T 


3 / -GTC(HEG)GACT 


(4t) 


38 


rGMP (5g) 


rC (RNA) 


3'-r(GUC(HEG)GACC) c 


(6c) 


14 



a Conditions: 0.5 mM hairpin in D 2 0 and 200 mM phosphate buffer, 400 mM NaCl, 80 mM MgC^, pH 8.9, uncorrected for deuterium effect, 20°C. Bold 
letters in hairpin sequences are templating bases; HEG, hexaethylene glycol linker. 
b Determined by fit. 

c 01igoribonucleotide (2 mM) and pH 7 to avoid hydrolysis of RNA under more basic conditions. 
d No saturation observed upon addition of up to 2000 eq. dCMP (1 M). 



extension reaction 



M+Pi ™ p 2 



(i) 



with M the monomer, Pi the primer-template complex, 
P2 the extended primer, and k the second order rate con- 
stant, is treated as a pseudo-first-order reaction with the ef- 
fective rate constant k! [M]. This approach is justified be- 
cause the monomer M was present in large excess. The 
unactivated nucleotide acting as inhibitor (Inh) and the 
primer-template complex are in a fast equilibrium with the 
inhibitor-primer-template complex Inh - Pi 



^Inh 

Inh + Pi ^ Inh 



(2) 



Therefore, the fraction of the free primer-template com- 
plex is reduced to 



[Pllfree 
[Pi] 



1 



1 + *i nh [Inh] 



(3) 



and the effective first order rate constant for the disappear- 
ance of Pi is given by 



/V[M] 



1 + K lnh [Inh] 1 + ^i nh [Inh] 



Taking the inverse, 

1 1 



^Inh 



[Inh] 



(4) 



(5) 



a linear relation between l//c e ff and the inhibitor concentra- 
tion is obtained. Plotting (1 //W) against [Inh] should yield a 
straight line with a slope Ki^/k and an intercept (l/k). Di- 
viding slope by intercept yields Ki n ^. For cases with a small 
number of data points and a strong inhibitory effect (at the 
highest concentration of the inhibitor), a large relative error 
in the intercept can result. Therefore, we used the actual ex- 
perimental values of k! [M] at [Inh] = 0 for such cases (see 
Supplementary Table SI and Supplementary Figures S23- 
S29 in the Supplementary material). 

Figure 5 shows the components of the kinetic assays 
performed. We used oligodeoxynucleotide sequences and a 
primer with a 3 / -terminal 3 / -amino-2 / ,3 / -dideoxynucleoside 
(40). Monomers were OAt-esters of 3 / -deoxynucleotides 



(35). A total of 16 different sequence motifs were employed, 
four each for A/A, C/C, G/G, or T/T as neighboring bases 
to the four templating bases (A, C, G or T). 

The rates of individual extension reactions were deter- 
mined in assays monitored by MALDI-ToF MS (57), using 
conditions that allow for quantitative detection of oligonu- 
cleotides (58). Figure 6 shows a representative data set for 
each of the four nucleotides (A, C, G and T). Additional 
data can be found in the Supplementary data (Supple- 
mentary Figures S15-S21). Primer extension was measured 
at 20° C in the absence or the presence of a downstream- 
binding 'helper' oligonucleotide that provides additional 
stacking interactions (35), whose effect may be similar 
to that of nucleotides in template-directed polymerization 
(59). Five additional data sets were acquired in the absence 
of the downstream-binding strand at 10°C under conditions 
known to give near-quantitative conversion (40). The dis- 
sociation constants for the complexes of the free deoxynu- 
cleotides and the oligonucleotide duplexes are compiled in 
Table 2. 

The extension rate for the monomer reacting the fastest 
(dGMP-OAt on the template with TCT as core motif) gave 
a t\ 12 of 2 min in the absence of any inhibitor. This value de- 
creased by a factor of 20 in the presence of 10 eq. of dGMP 
(Supplementary Figures S31 and 32, Supplementary ma- 
terial), indicating strong binding. For the extension with 
dAMP-OAt, directed by the TTT template motif, the ty 2 
increased from 7 to 60 min when adding 10 eq. of dAMP. 
There was also a significant effect of the neighboring nu- 
cleotides. When adenines were the neighbors of the incom- 
ing nucleotide (TNT templates, entries 1^ of Table 2), bind- 
ing was strongest. The absence of the downstream-binding 
strand led to an ~3- to 4-fold increase in at 20° C (entries 
l^r and 17-20, Table 2), whereas lowering the temperature 
from 20 to 10°C had a modest and less uniform effect on 
binding (see also Supplementary Figure S30 in the Supple- 
mentary data). 

Overall, the dissociation constants found range from 2 
mM for dGMP binding to a TCT template in the presence 
of a helper strand and 200 mM for TMP binding to a TAT 
template in the absence of a helper. They are thus close to 
the values found for hairpins by NMR (Table 1). The slow- 
est reaction (incorporation of T on an A-template in the ab- 
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Figure 5. Oligonucleotide sequences and nucleotides for template-directed primer extension reaction in the presence or absence of an unactivated (free) 
deoxynucleotide as inhibitor. Assays at increasing concentrations of inhibitor were performed in the presence or the absence of a downstream-binding 
oligonucleotide that provides additional stacking interactions to the incoming nucleotide. Conditions: 3.6 or 7.2 mM monomer, 0-72 mM free nucleotide, 
primer extension buffer (200 mM HEPES, 80 mM MgCl 2 , 400 mM NaCl, pH 8.9), 20°C. 
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Figure 7. Association constants for dNMPs binding to primer-template 
duplexes displaying complementary base at 20° C, as obtained by averaging 
over values of NMR titrations and inhibitory studies (first nine entries of 
Table 1 and entries 1-20 of Table 2). See Tables 1 and 2 for conditions. 
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Figure 6. Binding of deoxynucleotides to primer-template complexes re- 
veals itself through inhibition of primer extension. Kinetics of extension 
of primer 5 / -CGCACGA-NH 2 -3 / (8a) by OAt-esters of deoxynucleotides 
as templated by TNT-type sequences, where the templating base N is A, 
C, G or T at increasing concentrations of unactivated dNMP (la-t) as in- 
hibitor, in the absence of a downstream-binding oligonucleotide at 20° C. 
Conditions: 36 |jlM primer, 3.6 or 7.2 mM dNMP-OAt (7a-t), 200 mM 
HEPES buffer, pH 8.9, 400 mM NaCl, 80 mM MgCl 2 . Symbols are exper- 
imental data and lines are monoexponential fits. 



sence of a downstream-binding oligonucleotide) also gave 
the smallest inhibitory effect for the unactivated nucleotide 
(Figure 6a). Even in the presence of 10 eq. of unactivated 



TMP, only a modest slowdown of the extension was ob- 
served. In contrast, addition of 10 eq. of dNMP led to a very 
significant drop in rate for each of the remaining three bases 
(N = A, C or G; Figure 6b-d), confirming a significant oc- 
cupation of the extension site by the monomers. Here, un- 
like in the hairpin case, the binding of dAMP and dCMP 
was rather similar to that of dGMP (Figure 6 and Table 2). 
Overall, binding strength decreases in the following order 
G>A^C>>T. 

This order is also reflected in the global binding constants 
for the four different deoxynucleotides shown in Figure 7 
that were calculated by averaging over the data at 20° C for 
each given base (Tables 1 and 2). The purines bind most 
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Table 2. Dissociation constants (K^'s) for nucleotides binding to tem- 
plate:primer complexes, as determined by inhibitor kinetics from global 
fits to primer extension data at increasing concentrations of free nucleotide 
inhibitors. The letter in the middle of the template sequence denotes the 
templating base 



b) 



Entry 


Nucleotides a 
7a-t/la-t 


Template 
B^B 2 


8nnn 


dbo b 


Temp. 

(°Q 


(mM) 


i 
i 


T 
1 


TAT 


oiai 


_i_ 


_u 






r* 

V I 


TfT 
1L1 




_i_ 


90 


L 


J 


r 1 

V 


TPT 

1 VI 1 


Rfrrf 


_i_ 


1C\ 
_u 


C 


A 
*+ 


A 


TTT 
111 




_i_ 


90 


c 

J 


J 


T 
1 


vJ/vvJ 


8gag 


_i_ 


90 


OJ 


O 


VI 


vJv^vJ 


8gcg 


_i_ 


90 


1 0 
1U 


7 

/ 


V 


VjVJ Vj 


8ggg 


_i_ 


90 


1 8 


0 

o 


A 




sgtg 


_i_ 


90 


c 

J 


Q 


T 
1 


P AT 


8cac 


_i_ 


90 




1 n 


r" 

V I 


rrr 

V V. V 


8ccc 


_l_ 


90 


7 

/ 


i 1 
1 1 


r 1 

V 


rrr 

V^VJV^ 


8cgc 


_i_ 


90 


o 


12 


A 


CTC 


8ctc 


+ 


20 


7 


13 


T 


AAA 


8aaa 


+ 


20 


77 


14 


G 


ACA 


8aca 


+ 


20 


5 


15 


C 


AGA 


8aga 


+ 


20 


15 


16 


A 


ATA 


8ata 


+ 


20 


5 


17 


T 


TAT 


8tat 




20 


200 


18 


G 


TCT 


8tct 




20 


15 


19 


C 


TGT 


8tgt 




20 


19 


20 


A 


TTT 


8ttt 




20 


17 


21 


T 


CAC 


8cac 




10 


144 


22 


G 


CCC 


8ccc 




10 


5 


23 


C 


CGC 


8cgc 




10 


30 


24 


A 


CTC 


8ctc 




10 


7 


25 


C 


TGT 


8tgt 




10 


10 



a Monomer concentration: 3.6 mM for dAMP-OAt (7a), dCMP-OAt (7c) 
or dGMP-OAt (7g), and 7.2 mM for TMP-OAt (7t), and 0, 1, 5 or 10 eq. 
of unactivated monomer (la-t). 
b Downstream-binding oligonucleotide (dbo). 

c Dissociation constant for deoxynucleotide, calculated from K\ n h values. 



strongly, followed by dCMP, and TMP, which binds ~ 18- 
fold less strongly than guanine. This suggests that the com- 
bination of stacking interactions (strongest for the purines), 
number of hydrogen bonds and strength of secondary elec- 
trostatic interactions (60) governs binding strength, with the 
former providing a significant portion of the overall free en- 
ergy of binding. 



Binding isotherms 

With the binding constants in hand, we asked to what ex- 
tent the reaction site of primer extension was occupied by 
the cognate nucleotide at a given concentration. Figure 8 
shows calculated occupancies for different nucleotides and 
binding scenarios. When the reaction site is occupied, the 
primer terminus is protected from side reactions (38). An 
extension site occupied by the correctly paired monomer 
is also blocked from untemplated misincorporations. Un- 
templated reactions are common, and over-extension of 
primers, beyond the length of the template, are frequently 
observed (37,44). Finally, in the bound state, the activated 
nucleotide will be at least partially protected from side reac- 
tions with other nucleotides (unspecific polymerization, py- 
rophosphate formation, etc.) (61) and, being sterically less 
accessible, less prone to hydrolyze. 
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Figure 8. Occupancy of extension site by the deoxynucleotide comple- 
mentary to the templating base at 20° C, as calculated for different con- 
centrations of 2 / -deoxynucleotides la-t using binding constants reported 
in Table 1 or Table 2. Binding to (a) hairpins 2a, 2c, 2g or 4t, (b) template- 
primer duplexes 8tnt:9a-t and (c) template-primer duplexes 8tnt:9a-t in 
the presence of downstream-binding oligonucleotide lOa-t. Note the dif- 
ferent scales of the x-axes in (a) and (b)/(c). 



The binding isotherms shown in Figure 8 show that at 
100 mM nucleotide concentration, only dGMP achieves 
near-quantitative occupancy of the primer extension site. 
For TMP, reaching a similar occupancy level would re- 
quire unrealistically high concentrations. At the low mil- 
limolar concentrations typical for dNTPs in the cell, only 
basal binding occurs at the extension site in the absence of a 
helper oligonucleotide (or a polymerase). A neighboring nu- 
cleotide may provide a modest helper-like effect (59). Com- 
parison of Figure 8b and c shows how a strongly binding 
neighbor can help with the incorporation of a weakly pair- 
ing monomer. It is known that one weakly pairing base in 
a sequence can become a 'block' for enzyme-free copying 
(62). 

Binding and rate of extension 

Next, we asked how well binding correlates with rate. 
A strong correlation would suggest that the strength of 
the template effect is the dominant factor that determines 
whether enzyme-free primer extension occurs or not. Fig- 
ure 9 shows a plot of rates versus binding constants for the 
16 different template sequences studied (Table 2). It can be 
discerned that a loose correlation exists. Still the correla- 
tion is weak enough to suggest that other factors also play 
a role. Probably, the second step of the two-step mechanism 
proposed earlier for extension of amino-terminal primers 
(pseudorotation of a pentavalent intermediate; expulsion of 
the leaving group) (40) has a slightly different sequence de- 
pendence than the non-covalent binding equilibrium. 

Inhibition by spent monomers 

Since the binding constants for unactivated nucleotides 
are also inhibitory constants, the binding data provide a 
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Figure 9. Correlation between dissociation constants and rates of exten- 
sion (determined for 3.6 mM concentration of dAMP-OAt 7a, dCMP-OAt 
7c and dGMP-OAt 7g, or 7.2 mM concentration of TMP-OAt It in the 
presence of downstream-binding oligonucleotides lOa-t at 20°C) for the 
templating sequences listed in Table 2. The lines are linear functions ob- 
tained using regression analysis; solid black line, all 16 values (r 2 = 0.775), 
broken line: highest data point excluded (r 2 = 0.209). 



quantitative answer to the question of how important the 
inhibitory effect of spent monomers (nucleotides formed 
through hydrolysis of activate nucleotides) is for enzyme- 
free copying in different experimental scenarios. Figure 10 
shows representative kinetics of primer extension and hy- 
drolysis, together with the calculated occupancy of the ex- 
tension site, assuming that monomer and spent monomer 
bind with the same affinity. 

Figure 10a shows that for a highly reactive, amino- 
terminal primer and a monomer binding strongly, exten- 
sion is so fast that the reaction is complete before hydroly- 
sis produces a significant concentration of inhibitor. Figure 
10b shows the case of a highly reactive, amino-terminus but 
a weakly pairing monomer (TMP-OAt 7t binding to hair- 
pin 2a). Here, the extension is so slow that the formation 
of a significant concentration of spent monomer occurs be- 
fore the reaction is over. But, the extent of binding of the 
inhibitor is minimal, so that its formation is inconsequen- 
tial. Finally, Figure 10c shows the case of a less reactive 
RNA primer, combined with strongly binding GMP-OAt 
as monomer. Here, extension is so slow that hydrolysis can 
catch up with the desired copying process, so that the in- 
hibitor formed does block extension significantly. (Perhaps, 
such an inhibitory effect can help to smooth out the differ- 
ences in reaction rates between strongly and weakly pair- 
ing nucleotides, making it more likely that the weakly bind- 
ing nucleotides compete successfully with the more strongly 
pairing ones.) 
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Figure 10. The role of inhibition in different chemical primer extension 
scenarios. On the right of each part, kinetics of primer extension (black 
dots) and monomer hydrolysis (open circles) are shown, and on the left, 
the calculated occupancy of the extension site by monomer or hydrolyzed 
monomer (inhibitor) at 0 min is shown as a bar graph, (a) Amino-terminal 
primer 9a (36 |jlM), template 8tct, helper strand 10a and 3.6 mM dGMP- 
OAt (7g); (b) 3 / -amino-terminal hairpin 2a (36 |jlM) and 3.6 mM dTMP- 
OAt (It); (c) RNA primer 13g (100 |jlM), template 12ccg, helper strand 14c 
and 20 mM rGMP-OAt (15g). The values used for calculating occupan- 
cies are 2 mM for dGMP, 260 mM for dTMP and 14 mM for GMP (Tables 
1 and 2). See also the Supplementary Material (Supplementary Figures 
SI 5b, S22 and S33). Kinetics and hydrolysis data for the RNA case are 
from reference (29). 



Yield of primer extension 

Incomplete conversion is a key obstacle to enzyme-free 
replication of oligonucleotides (26). Therefore, we asked to 
what extent the yield of enzyme-free copying (and thus the 
perhaps most critical step of spontaneous replication, other 
than strand separation) can now be predicted, based on the 
binding constants and rates for the extension and the hy- 
drolysis of monomers. First, we focused on an RNA-based 
system that is mechanistically simpler than those involving 
amino-terminal primers (40). All-RNA systems are at the 
focus of studies on prebiotic evolution (16,63). Figure 11 
shows the model used for the simulations. 



M + Pi M..P! — ► M-P-, — ► P 2 



lnh + P 1 ^=^ Inh.Pi 
*d 

Figure 11. Binding equilibria and reactions for extension of a primer (Pi) 
by an activated nucleotide as monomer (M). Hydrolysis of the monomer 
produces a free nucleotide that acts as an inhibitor (Inh). The non-covalent 
binding of both activated and free nucleotide to the primer are governed by 
the dissociation constant (K^), while the rates of extension (k C0Y ) and hy- 
drolysis (&h) govern the fate of the monomer. It is assumed that the primer 
is stably bound to the template. 
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Figure 12. Oligonucleotide sequences and ribonucleotides for template- 
directed RNA primer extension reaction in the presence or absence of an 
unactivated (free) ribonucleotide as inhibitor, as described (29). 



As shown in detail in the Supplementary material, the 
yield of product P2 can be calculated by Equation (6): 



Y jn = 1 — exp 



k cov [M] 0 (Qxp[-k h t] - 1) 
L k h (K d + [M] 0 + [Inh]) . 



(6) 



Yields of extension of RNA primers were calculated us- 
ing the binding constant determined for GMP and RNA 
hairpin 5g (14 mM, Table 1) and the data for extension 
and hydrolysis recently reported by Deck et al for the se- 
quence system shown in Figure 12 (29). Assuming that ac- 
tivated and unactivated nucleotide have the same affinity 
for the template and that the monomer binding equilibrium 
is reached rapidly, the occupation number is given by a = 
[M]o/(^d + [M]o)- Using the experimental K d value of 14 
mM (Table 1), a = 0.59 for assays performed at 20 mM 
monomer concentration. 

At early time points, the kinetics are not yet affected by 
hydrolysis of the monomer. In this case, the effective first or- 
der rate of the reaction can be determined using k = ak COY . 
A mono-exponential fit to the first four experimental data 
points shown in Figure 13a yields k = 0.16 h _1 leading 
to k cow = 0.27 h _1 . This value is similar to the rate of the 
extension with GMP-OAt on a slightly different template 
motif (UCU, where C is the templating base), which occurs 
with 0.38 h _1 under the same experimental conditions (29). 
Either value is close to the rate of hydrolysis under assay 
conditions (k^ = 0.15 h _1 ) (29), as expected for the 'third 
scenario' of Figure 10c. 

Figure 13a shows plots of the simulated progress of 
primer extension during the first two days reaction time 
of the RNA extension and four different regimes. The first 
regime is that of a hypothetical extension without any in- 
hibitory effect of hydrolyzed monomer (dashed black line). 
The second (solid black line) was calculated with the full 
model that takes concomitant formation of the inhibitor 
through hydrolysis into account. The third and the fourth 
case (dark gray and pale gray line, respectively) were calcu- 
lated for extensions in the presence of 5 or 20 mM GMP 
(5g) (29). 

For extension without addition of inhibitor, the calcu- 
lated curves are close to the available experimental data 
(black circles in Figure 13), but the drop in yield induced 
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Figure 13. Simulated time-dependent yield of extended RNA primer (Fig- 
ure 12) with and without hydrolysis and inhibition at initial concentra- 
tions of 20 mM monomer (GMP-OAt) and 268 |jlM primer, (a) Symbols 
represent published experimental values (29) for assays with different con- 
centrations of inhibitor (GMP) added; circles: none, triangles 5 mM, dia- 
monds 20 mM GMP. Solid lines are calculated using Equation (6). Values 
used for the simulation: K d = 14 mM; k cov = 0.27 h -1 , k h = 0.147 h -1 . 
The dashed black line shows hypothetical kinetics without formation of in- 
hibitor through hydrolysis, (b) Same as (a) but with a value increased 
by a factor of 3, as expected for this system containing a downstream- 
binding oligonucleotide. 



by addition of the inhibitor is represented in a qualitative 
sense only. For the simulation of Figure 13a, the binding 
constant from the hairpin system (Figure 2) was used, even 
though the experimental data was from a full extension sys- 
tem with downstream-binding helper oligonucleotide (Fig- 
ure 12). Helper oligonucleotides typically give a 3-fold de- 
crease in Kd, at least in the better studied DNA case (com- 
pare the first four entries of Table 2 with those in the lower 
part of the table where the same sequence motifs were as- 
sayed without helper). When the unactivated GMP was as- 
signed a 3-fold lower value than the activated monomer, 
a near-perfect agreement of calculated and experimental 
data was obtained (Figure 13b). 

Finally, for assays with amino-terminal primers and 
DNA templates, the mechanistic situation is more compli- 
cated, and inhibition does not play a significant role. To test 
the limits of our theory, we simulated assays with decreas- 
ing concentration of the monomer, down to the micromo- 
lar range (10 eq. or even 1 eq. of activated monomer). This 
case is more challenging, as the k cov values are more diffi- 
cult to extract for biphasic kinetics, and because at so low 
a monomer concentration, side reactions from trace impu- 
rities, such as residual acetate, become more prevalent (38). 
Figure 14 shows the results for A and G as monomers. 

It can be discerned that Equation (6) predicts the time- 
and concentration-dependent yields well for all but the very 
lowest concentrations of the monomers. Figure 14c shows 
that assuming a loss of reactive species, so that the effective 
concentration of the monomer is half of what is assumed 
in the ideal case, suffices to get a satisfactory agreement be- 
tween theoretical and experimental data, even at the very 
lowest monomer concentration. Apparently, binding con- 
stants, global rate constants for the covalent step(s), and 
rate constants for hydrolysis of monomers largely suffice to 
explain incomplete conversion, even in this case. 
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Figure 14. Simulated time-dependent yield of extension of amino- 
terminal DNA primer 9g at decreasing concentration of monomer in the 
presence of downstream-binding oligonucleotide lOg. Lines are calculated 
using Equation (6) with k cov values obtained from rate constants (assuming 
monoexpoential kinetics) and occupation numbers based on dissociation 
constants for nucleotides. Filled circles are experimental data for monomer 
concentrations of 3.6 mM (red), 0.36 mM (blue), 0.18 mM (green) or 0.036 
mM (pink), (a) Monomer is A (dAMP-OAt) and template is 8ctc; values of 
the simulation: K d = 6.9 mM; k/ t = 0. 109 h _1 , k cov = 8.3 h _1 . (b) Monomer 
is G (dGMP-OAt) and template is 8ccc; with values of = 6.8 mM; kh = 
0.093 h _1 and k cov = 12.2 h _1 . (c) Same as (b), except that a 2-fold lower 
monomer concentration was assumed for the pink data at the end of the 
dilution series (0.018 mM). 



DISCUSSION 

Our results show that the dissociation constants for 
nucleotide-primer/template complexes are in the millimo- 
lar range. They are generally weaker than previously 
thought (51). The sequence dependence shows similarities 
to that found in the study on the rates of chemical primer ex- 
tension (40). A large difference in binding strength is found 
between thymine and the other three bases, amounting to 
approximately one order of magnitude difference in disso- 
ciation or binding constants. Thymidine monophosphate is 
the only base for which near- saturation of the extension 
site is unrealistic at room temperature, even at concentra- 
tions approaching the physical limit. Lowering the temper- 
ature to 10°C does not eliminate this problem, nor does the 
presence of a downstream-binding oligonucleotide, though 
it does have a significant effect (Figure 15). Lowering the 
temperature further can be expected to tighten binding, and 
thus to improve yields (28). 

Even in the best of all templating environments, with 
two large, hydrophobic adenine bases as neighbors (at the 
primer terminus and the 5 / -end of the helper strand), the K& 
value for TMP remains as high as 59 mM (first entry, Table 
2). Having purines as stacking partners at the neighboring 
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Figure 15. Effect of nucleobase and sequence context on binding of nu- 
cleotides to templating bases: heat map representation of representative 
dissociation constants for complexes between deoxynucleotides and ter- 
mini of hairpins, primer-template complexes or primer-template com- 
plexes with downstream-binding oligonucleotide at 20° C. The color bar 
on the right-hand side is a graphical definition of how color intensity codes 
for binding strength. Data are from Tables 1 and 2. 



positions helps all incoming deoxynucleotides, but the over- 
all effect of binding in a TNT' template region or binding 
to a less favorable 4 GNG' sequence, with two cytosines as 
stacking partners, changes the stability of the complex by 
a factor of five at best, and as little as a factor of 1.4 in the 
case of dAMP. Being without any downstream interactions, 
as in the case of binding a base at the very terminus of a 
template/hairpin, leads to a decrease in binding that can be 
as large of 7-fold, as in the case of dCMP binding to hairpin 
2g (Table 1) when compared to the long template 8tgt in the 
presence of a helper strand (Table 2). For the most tightly 
binding deoxynucleotide, dGMP, the dissociation constants 
found range from 2 mM on template 8tct with helper and 
16 mM for hairpin 3c. (This concentration range is close to 
that of the intracellular concentration of nucleoside triphos- 
phates (NTPs) found in present-day cells (64).) 

The results of the exploratory experiments on mis- 
matched binding provide a glimpse of the energetic conse- 
quence of mismatches. A G:T wobble base pair, most prob- 
ably the most problematic of all mismatches (42,43) leads 
to binding that is just 10-fold weaker than that of the corre- 
sponding G:C combination (templates 4t and 2c, Table 1). 
Its K& value is lower than that of the canonical T:A pairing 
(lt:2a), emphasizing how difficult it is to suppress this bind- 
ing mode with natural T as templating base (65) at equimo- 
lar concentrations of all four monomers. On the other hand, 
no binding was detectable in the NMR titration for mis- 
matched deoxynucleotide lc and hairpin 4t, demonstrating 
that not all mismatches are problematic. Also, the pairing 
between GMP (5g) and hairpin 6c consisting of RNA gave 
a value of 14 mM, which is close to the values found 
for DNA hairpins 2c and 3c. This suggests that backbone 
structure has a minor effect on the strength of base pairing 
in this sequence context. 
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The nature of the primer (amino-terminal DNA versus 
RNA) does make a significant difference in terms of the 
importance of inhibition, though. As Figures 10 and 13 
show, extension of an RNA primer, with slow-reacting nat- 
ural backbone, suffers significantly from competitive inhi- 
bition by unactivated monomer, whereas the more reactive 
amino-terminal primer forming phosphoramidate linkages 
does not. This insight had eluded us and others, when ki- 
netic data was available only. 

Figure 8 allows one to gauge whether further increases 
in monomer concentration can be expected to improve 
yields for a given setting or whether this will be futile be- 
cause near-saturation of the extension site has already been 
achieved. The data also helps to understand why submil- 
limolar monomer concentrations require re-activation of 
spent monomers to achieve high yields in chemical primer 
extension (38). 

Finally, the present data also allows a first glimpse at how 
the active site of polymerases improves binding of nucleo- 
side phosphates. For example, the complex of dGTP and 
the phi29 DNAP polymerase was recently reported to form 
with a K& of 1 .4 (xM (66), a value that is just three orders of 
magnitude smaller than some of the K&s measured here for 
dGMP (and part of that increase in affinity is probably due 
to the interactions between the additional pyrophosphate 
and the Mg 2+ ions in the active site). 



CONCLUSIONS 

Our manuscript reports methodologies for measuring bind- 
ing of nucleotides to templating bases, binding data and 
a model for calculating yields of chemical primer exten- 
sions. The model has been validated by simulating incom- 
plete copying reactions for which experimental data are 
available. As Figures 13 and 14 show, we have also success- 
fully separated binding from intrinsic reactivity in the ac- 
tive extension complex (the k cov values show that amino- 
terminal primers are ~ 3 0-fold more reactive than RNA 
primers). Our approach should also allow the quantita- 
tive prediction of the inhibitory effect of unactivated or 
spent monomers on enzyme-free primer extension for other 
chemistries and assay conditions, provided that binding 
constants and kinetic constants for extension and monomer 
hydrolysis are available. All three types of constants are ac- 
cessible by straightforward experiments, using NMR (bind- 
ing constant and hydrolysis) and extension assays moni- 
tored by gel electrophoresis or mass spectrometry. 

Significant issues remain before enzyme-free replication 
of nucleic acids may be shown experimentally (67). Quanti- 
tative simulations of replication scenarios probably should 
include a matrix of fidelity values and quantitative data on 
the stalling after misincorporations (41). Further, it is in- 
teresting to ask how scenarios involving several monomers 
pairing simultaneously with a longer stretch of template 
affect reactivity, or how combinations of monomer exten- 
sion and ligation, compared to the purely monomer-based 
regime known from present-day replication and transcrip- 
tion, perform. Efforts to tackle such systems are under way 
in our laboratories. 



SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online, includ- 
ing protocols, NMR data, primer extension data, a more 
elaborate treatment of inhibitor kinetics, kinetics of the hy- 
drolysis of monomers, and description of the model for pre- 
dicting yields. 

ACKNOWLEDGMENTS 

The authors thank Annette Hochgesand for expert tech- 
nical assistance during an early phase of this project, Dr 
Christopher Deck for sharing data and Jan Meisner for per- 
forming early computational work and helpful discussions. 

FUNDING 

Deutsche Forschungsgemeinschaft [RI 1063/8-2 to C.R.]. 
Funding for open access charge: Deutsche Forschungsge- 
meinschaft (DFG) . 

Conflict of interest statement. None declared. 
REFERENCES 

1. Kornberg,A. and Baker,T.A. (2005) DNA Replication, 2nd edn. 
University Science Books, Mill Valley, CA. 

2. Woese,C.R. (1967) The Genetic Code; The Molecular Basis For 
Genetic Expression. Harper & Row, New York, NY. 

3. Kunkel,T.A. and Bebenek,R. (2000), DNA replication fidelity. Annu. 
Rev. Biochem., 69, 497-529. 

4. Bryant,F.R., Johnson,K.A. and Benkovic,S.J. (1983) Elementary 
steps in the DNA-polymerase-I reaction pathway. Biochemistry, 22, 
3537-3546. 

5. Doublie,S., Tabor,S., Long,A.M., Richardson,C.C. and 
Ellenberger,T. (1998) Crystal structure of a bacteriophage T7 DNA 
replication complex at 2.2 A resolution. Nature, 391, 251-258. 

6. Lieberman,K.R., DahlJ.M., Mai,A.H., Cox,A., Akeson,M. and 
Wang,H. (2013) Kinetic mechanism of translocation and dNTP 
binding in individual DNA polymerase complexes. J. Am. Chem. 
Soc, 135, 9149-9155. 

7. PiccirilliJ.A., Krauch,T., Moroney,S.E. and Benner,S.A. (1990) 
Enzymatic incorporation of a new base pair into DNA and RNA 
extends the genetic alphabet. Nature, 343, 33-37. 

8. Kool,E.T. (2002) Replacing the nucleohases in DNA with designer 
molecules. Acc. Chem. Res., 35, 936-943. 

9. Krueger,A.T. and Kool,E.T. (2007) Model systems for understanding 
DNA base pairing. Curr. Opin. Chem. Biol., 11, 588-594. 

10. Malyshev,D.A., Pfaff,D.A., Ippoliti,S.L, Hwang,G.T., Dwyer,T.J. and 
Romesberg,F.E. (2010) Solution structure, mechanism of replication, 
and optimization of an unnatural base pair. Chem. Eur. J., 16, 
12650-12659. 

11. Moran,S., Ren,R.X.E, Rumney,S. and Kool,E.T. (1997) 
Difluorotoluene, a nonpolar isostere for thymine, codes specifically 
and efficiently for adenine in DNA replication. J. Am. Chem. Soc. , 
119, 2056-2057. 

12. Egli,M. (2012) The steric hypothesis for DNA replication and 
fluorine hydrogen bonding revisited in light of structural data. Acc. 
Chem. Res., 45, 1237-1246. 

13. Hwang,G.T. and Romesberg,F.E. (2006) Substituent effects on the 
pairing and polymerase recognition of simple unnatural base pairs. 
Nucleic Acids Res., 34, 2037-2045. 

14. Orgel,L.E. and Lohrmann,R. (1974) Prebiotic chemistry and 
nucleic-acid replication. Acc. Chem. Res., 7, 368-377. 

15. Kaiser, A. and Richert,C. (2013) Nucleotide-based copying of nucleic 
acid sequences without enzymes. J. Org. Chem., 78, 793-799. 

16. Orgel,L.E. (2004) Prebiotic chemistry and the origin of the RNA 
world. Crit. Rev. Biochem. Mol. Biol, 39, 99-123. 

17. AtkinsJ.E, Gesteland,R.F. and Cech,T. (eds.). (2011) RNA Worlds. 
Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY. 



Nucleic Acids Research, 2014, Vol 42, No. 11 7419 



18. Shechner,D.M., Grant,R.A., Bagby,S.C, Koldobskaya,Y, 
Piccirilli,!A. and Bartel,D.R (2009) Crystal structure of the catalytic 
core of an RNA-polymerase ribozyme. Science, 326, 1271-1275. 

19. Wochner,A., Attwater,!, Coulson,A. and Holliger,P. (2011) 
Ribozyme-catalyzed transcription of an active ribozyme. Science, 
332,209-212. 

20. Sulston,!, Lohrmann,R., Orgel,L.E. and Miles,H.T. (1968) 
Nonenzymatic synthesis of oligoadenylates on a polyuridylic acid 
template. Proc. Natl. Acad. Sci. U.S.A., 59, 726-733. 

21. Weimann,B.J., Lohrmann,R., Orgel,L.E., Schneider-Bernloehr,H. 
and Sulston, IE. (1968) Template-directed synthesis with 
adenosine-5'-phosphorimidazolide. Science, 161, 387. 

22. Lohrmann,R., Bridson,P.K. and Orgel,L.E. (1980) Efficient 
metal-ion catalyzed template-directed oligonucleotide synthesis. 
Science, 208, 1464-1465. 

23. Inoue,T. and Orgel,L.E. (1983) A nonenzymatic RNA polymerase 
model. Science, 219, 859-862. 

24. Kurz,M., G6bel,K., Hartel,C. and G6bel,M.W. (1997) Nonenzymatic 
oligomerization of ribonucleotides on guanosine-rich templates: 
suppression of the self-pairing of guanosine. Angew. Chem. Int. Ed, 
36, 842-845. 

25. Hill,A.R., Orgel,L.E. and Wu,T.F. (1993) The limits of 
template-directed synthesis with 

nucleoside-5'-phosphoro(2-methyl)imidazolides. Orig. Life Evol. 
Biospk, 23, 285-290. 

26. Zielinski,M., Kozlov,I.A. and Orgel,L.E. (2000) A comparison of 
RNA with DNA in template-directed synthesis. Helv. Chim. Acta, 83, 
1678-1684. 

27. Vogel,S.R., Deck,C. and Richert,C. (2005) Accelerating chemical 
replication steps of RNA involving activated ribonucleotides and 
downstream-binding elements. Chem. Commun. , 2005, 4922-4924. 

28. Vogel,S.R. and Richert,C. (2007) Adenosine residues in the template 
do not block spontaneous replication steps of RNA. Chem. 
Commun.,2007, 1896-1898. 

29. Deck,C, Jauker,M. and Richert,C. (2011) Efficient enzyme-free 
copying of all four nucleobases templated by immobilized RNA. Nat. 
Chem.,?,, 603-608. 

30. Adamala,K. and Szostak,! W. (2013) Nonenzymatic 
template-directed RNA synthesis inside model protocells. Science, 
342, 1098-1100. 

31. Lehman,N. and Hayden,E.J. (2011) Template-directed RNA 
polymerization: the taming of the milieu. ChemBioChem, 12, 
2727-2728. 

32. Szostak,! W. (2012) The eightfold path to non-enzymatic RNA 
replication. J. Syst. Chem. , 3, 2. 

33. Lohrmann,R. and Orgel,L.E. (1976) Template-directed synthesis of 
high molecular weight polynucleotide analogues. Nature, 261, 
342-344. 

34. Zielinski,W.S. and Orgel,L.E. (1985) Oligomerization of activated 
derivatives of 3'-amino-3'-deoxyguanosine on poly(C) and poly(dC) 
templates. Nucleic Acids Res., 13, 2469-2484. 

35. Hagenbuch,R, Kervio,E., Hochgesand,A., Plutowski,U. and 
Richert,C. (2005) Chemical primer extension: efficiently determining 
single nucleotides in DNA. Angew. Chem. Int. Ed, 40, 6588-6592. 

36. R6thlingsh6fer,M., Kervio,E., Lommel,T, Plutowski,U., 
Hochgesand,A. and Richert,C. (2008) Chemical primer extension in 
seconds. Angew. Chem. Int. Ed, 47, 6065-6068. 

37. Mansy,S.S., SchrumJ.R, Krishnamurthy,M., Tobe,S., Treco,D.A. 
and Szostak, J.W. (2008) Template-directed synthesis of a genetic 
polymer in a model protocell. Nature, 454, 122-125. 

38. R6thlingshdfer,M. and Richert,C. (2010) Chemical primer extension 
at submillimolar concentration of deoxynucleotides. J. Org. Chem., 
75, 3945-3952. 

39. Kaiser,A., Spies, S., Lommel,T. and Richert,C. (2012) Chemical 
primer extension at submillimolar concentration of deoxynucleotides. 
Angew. Chem. Int. Ed, 51, 8299-8303. 

40. Kervio,E., Hochgesand,A., Steiner,U.E. and Richert,C. (2010) 
Templating efficiency of naked DNA. Proc. Natl. Acad. Sci. U.S.A., 
107, 12074-12079. 

41. Rajamani,S., IchidaJ.K., Antal,T, Treco,D.A., Leu,K., 
Nowak,M.A., Szostak,! W. and Chen,I.A. (2010) Effect of stalling 
after mismatches on the error catastrophe in nonenzymatic nucleic 
acid replication. J. Am. Chem. Soc, 132, 5880-5885. 



42. Leu,K., Obermayer,B., Rajamani,S., Gerland,U. and Chen,I.A. 
(201 1) Fast and accurate nonenzymatic copying of an RNA-like 
synthetic genetic polymer. Nucleic Acids Res., 39, 8135-8147. 

43. Leu,K., Kervio,E., Obermayer,B., Turk-MacLeod,R.M., Yuan,C, 
Luevano,J.-M. Jr, Chen,E., Gerland,U., Richert,C. and Chen,I.A. 
(2013) Cascade of reduced speed and accuracy after errors in 
enzyme-free copying of nucleic acid sequences. J. Am. Chem. Soc, 
135, 354-366. 

44. SchrumJ.R, Ricardo,A., Krishnamurthy,M., BlainJ.C. and 
Szostak,! W. (2009) Efficient and rapid template-directed nucleic acid 
copying using 2'-amino-2', 

3'-dideoxyribonucleoside-5'-phosphorimidazolide monomers. J. Am. 
Chem. Soc, 131, 14560-14570. 

45. Griesang,N, Giessler,K., Lommel,T. and Richert,C. (2006) Four 
color, enzyme-free interrogation of DNA sequences with chemically 
activated, 3'-fluorophore-labeled nucleotides. Angew. Chem. Int. Ed, 
45, 6144-6148. 

46. GieBler,K., Griesser,H., G6hringer,D., Sabirov,T. and Richert,C. 
(2010) Synthesis of 3'-BODIPY-labeled active esters of nucleotides 
and a chemical primer extension assay on beads. Eur. J. Org. Chem., 
19, 3611-3620. 

47. Kyogoku,Y, Lord,R.C. and Rich,A. (1969) An infrared study of 
hydrogen-bonding specificity of hypoxanthine and other nucleic acid 
derivatives. Biochim. Biophys. Acta, 179, 10-17. 

48. Weisz,K., Jaehnchen,! and Limbach,H.-H. (1997) Low-temperature 
NMR studies on the structure of uridine dimers. J. Am. Chem. Soc, 
119, 6436-6437. 

49. Guerra,C.F, Bickelhaupt,F.M., SnijdersJ.G. and Baerends,E.J. 
(2000) Hydrogen bonding in DNA base pairs: reconciliation of 
theory and experiment. J. Am. Chem. Soc, 122, 4117-4128. 

50. Fonseca Guerra,C. and Bickelhaupt,F.M. (2002) Orbital interactions 
in strong and weak hydrogen bonds are essential for DNA 
replication. Angew. Chem. Int. Ed., 41, 2092-2095. 

51. Poater,!, Sart,M., Fonseca Guerra,C. and Bickelhaupt,F.M. (2011) 
Selectivity in DNA replication. Interplay of steric shape, hydrogen 
bonds, pi-stacking and solvent effects. Chem. Commun., 2011, 
7326-7328. 

52. Eisenhuth,R. and Richert,C. (2009) Convenient syntheses of 
3'-amino-2', 3'-dideoxynucleosides, their 5' -monophosphates, and 
3'-aminoterminal oligodeoxynucleotide primers. J. Org. Chem., 74, 
26-37. 

53. Altmann,S., Labhardt,A.M., Bur,D., Lehmann,C, Bannwarth,W., 
Billeter,M., Wuthrich,K. and Leupin,W. (1995) NMR studies of 
DNA duplexes singly cross-linked by different synthetic linkers. 
Nucleic Acids Res., 23, 4827-4835. 

54. Wijmenga,S.S., Mooren,M.M.W. and Hilbers,C.W. (1993) NMR of 
nucleic acids; from spectrum to structure. In: Roberts, GCK (ed), 
NMR of Macromolecules — A Practical Approach. IRL Press, 
Oxford, pp. 217-288. 

55. Fielding,L. (2000) Determination of association constants (Ka) from 
solution NMR data. Tetrahedron, 56, 6151-6170. 

56. Hirose,K. (2001) A practical guide for the determination of binding 
constants. J. Inch Phenom. Macrocycl. Chem., 39, 193-209. 

57. Rojas StutzJ.A. and Richert,C. (2001) A steroid cap induces 
enhanced selectivity and rates in non-enzymatic single nucleotide 
extensions of an oligonucleotide. J. Am. Chem. Soc. , 123, 
12718-12719. 

58. Sarracino,D. and Richert,C. (1996) Quantitative MALDI-TOF 
spectrometry of oligonucleotides and a nuclease assay. Bioorg. Med. 
Chem. Lett., 6, 2543-2548. 

59. Kanavarioti,A., Bernasconi,C.F. and Baird,E.E. (1998) Effects of 
monomer and template concentration on the kinetics of 
nonenzymatic template-directed oligoguanylate synthesis. J. Am. 
Chem. Soc, 120, 8575-8581. 

60. Pranata,!, Wierschke,S.G. and Jorgensen,W.L. (1991) OPLS 
potential functions for nucleotide bases-relative association constants 
of hydrogen-bonded base pairs in chloroform. J. Am. Chem. Soc. , 
113, 2810-2819. 

61. Hill,A.R. Jr, Wu,T. and Orgel,L.E. (1993) The limits of 
template-directed synthesis with 

nucleoside-5 / -phosphoro(2-methyl)imidazolides. Orig. Life Evol. 
Biospk, 23, 285-290. 



7420 Nucleic Acids Research, 2014, Vol 42, No. 11 



62. Rojas StiitzJ.A., Kervio,E., Deck,C. and Richert,C. (2007) Chemical 
primer extension-individual steps of spontaneous replication. Chem. 
Biodiv.,4, 784-802. 

63. Gesteland,R.R, Cech,T.R. and Atkins, J.F. (eds.). (2006) The RNA 
World, 3rd edn. Cold Spring Harbor Laboratory Press, Cold Spring 
Harbor, NY. 

64. Traut,T.W. (1994) Physiological concentrations of purines and 
pyrimidines. Mol. Cell. Biochem., 140, 1-22. 

65. Zhang, S., BlainJ.C, Zielinska,D., Gryaznov,S.M. and SzostakJ.W. 
(2013) Fast and accurate nonenzymatic copying of an RNA-like 



synthetic genetic polymer. Proc. Natl. Acad. Sci. U.S.A., 110, 
17732-17737. 

66. Lieberman,K.R., DahlJ.M., Mai,A.H., Cox,A., Akeson,M. and 
Wang,H. (2013) Kinetic mechanism of translocation and dNTP 
binding in individual DNA polymerase complexes. J. Am. Chem. 
Soc, 135, 9149-9155. 

67. Service,R.F. (2013) The life force. Science, 342, 1032-1034. 



